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Abstract 

We discuss the propagation of laser radiation through a medium of quantum prepared A-type 
atoms in order to enhance the insight into the physics of QSPT generator suggested in Phys. Rev. 
A 80, 035801 (2009). We obtain analytical results which give a qualitatively corerct description of 
the outcoming series of ultrashort optical pulses and show that for the case of alkali vapor medium 
QS-PT generation may be implemented under ordinary experimental conditions. 
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I. INTRODUCTION 



Atomic coherence is a widely accepted resource in atomic and laser physics for modifying 
the response of coherently prepared atomic systems and thereby canceling or enhancing the 
light absorption/reflection near atomic resonance [1-4]. Quantum prepared medium in fact is 
essentially a new state of matter, named in [5] phaseonium. The basic idea behind almost all 
such models is to modify the emissive and absorptive profiles and make them non-reciprocal. 
Among most important applications of quantum prepared medium, both in theory and 
experimental science,one should list adiabatic population transfer [6-9], amplification and 
lasing without inversion [10-13], sub-recoil laser cooling [14-17], light storage and stopping 
[18-21], atomic spectroscopy [22-26], etc. 

On the other hand, many problems in atomic physics, such as precision optical frequency 
metrology, excitation spectroscopy and quantum state engineering [27-30] , require a train of 
high-energy and high-repetition optical pulses. The most convenient method to accomplish 
this requirement today is the acousto-optic modulator [31], but still new principles and 
schemes are suggested in literature. A promising approach to this end is the usage of 
quantum coherency created between atomic levels in both bare and dressed approximations. 
For example, a considerable attention has been drawn to the parametric beating of a weak 
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probe field in the scheme of stimulated Raman scattering [32,33]. In this paper we suggest 
a new step in this direction. 

We consider propagation of an incident wave through a medium of A-type three-level 
atoms with a preliminary superposition between the low lying doublet states, a simplest 
construction for experimental realizations. We carry a detailed analysis of how the super- 
position nature of atomic states splits the incident monochromatic wave into a train of 
ultrashort (and high power) pulses in a three-level medium. Such a qualitatively new and 
dramatic effect of superposition principle has already been identified in our previous papers 
[34,35]. However, it has proceeded from a two-level medium and two (probe and pump) 
optical fields, thus somewhat disguising the inherent ability of superposition principle in 
this problem. Simultaneously, the technical realization of originally proposed scheme seems 
stringent as much as it is difficult to onset the preliminary quantum superposition between 
the energy levels in the optical range of Bohr frequencies. Presented in this paper scheme 
clarifies the main idea of [34,35] and extends its formalism into much headmost ideological 
and experimental conditions using A-type interaction scheme and only one travelling wave. 



II. ADIABATIC STATES OF A-TYPE ATOM 



The basic setup of interaction used throughout this article is presented in Fig. 1. The 
lower atomic states 1 and 2 are separated in energy by E a2 — E ai = hA , and are coupled 
with the upper state by a plane optical field of frequency u: 



E* = i^o(<2) t) exp(iujt — ikz) + ~E*q(z, t) exp(—icot + ikz), 



(1) 



where k = oo/c is the vacuum wave number. 

The amplitude E*o(z,t) is assumed to vary slowly compared to the rapidly oscillating ex- 



ponential functions, i.e. 



dE* (z,t)/dt « w~3 (z,t) and d~3 (z,t)/dz « k~E* (z,t) 
(Later we will take it as a constant when entering the gas medium, which is the most con- 
venient form for implication purposes). We also assume that the necessary conditions to 
preserve the coherency are ensured. Dephasing from the upper state is eliminated by assum- 
ing a far off-resonant interaction, while the relaxation or dephasing between ground levels 
1 and 2 lasts milliseconds and is sufficient for the objectives of this paper. Thus we will 
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FIG. 1: Energy-level diagram for the resonant single-photon process. States |al) and \a2) are of 
like parity whereas the upper state \b) is of opposite parity. Energy distance hAo between the 
doublet states is much smaller than the photon energy hu>. 



discuss the atom-radiation coupling staying in the limits of Schrodinger equation 



(2) 



The interaction is taken in electric dipole approximation, H and d are free atom Hamilto- 
nian and dipole moment operator respectively, and the atomic motion isn't included in the 
discussion. 

We look for the solution of Eq.2 in a form 



t) = ai(t)^ x (T) exp(-iE ai t/h)+a 2 (t)il> 2 ( t ) exp(-iE a2 t/h) + b(t)<p(T) exp(-iE b t/h), 



tively. E ai 2 and Eb are corresponding energies. After substituting ( |2j) into (jSJ) we 
apply the RWA neglecting terms with rapidly oscillating factor exp(±i(u + Ui j2 )t) with 
u)i j2 = Eb — E ai2 /h. We want to describe the phenomena under conditions where the reac- 
tive part of the atom-field coupling is predominant over the dissipative part so we proceed 
working in the range of resonance detunings much larger than the Doppler and homogeneous 
widths of optical transitions. This simultaneously gives the possibility of first order adia- 
batic elimination procedure, algebraically presenting the excited state amplitude through 
the lower lying ones: 




b(t) 




exp(— i{oj— uj\)t)a\{t) 




exj)(-i(u—U2)t)a 2 {t). (4) 



h{uj — oji) 



h(ui - oj 2 ) 



Here ~~ct i and d " 2 are al — 6 and a2 — 6 optical transition matrix elements correspondingly. 
For simplicity of later calculations we redefine the level 2 probability amplitude 

A 2 (t) = a 2 (t)exp(-zA t), A x {t) = a x {t), (5) 

and rescale the problem parameters to dimensionless ones: 

T = A t, — =UJ-U h — =U-W2, (6) 

^0 ^0 

, f \ rfi^o(^r) ~t '2^0 (z,t) 

The doublet amplitudes now will take a relatively simple form 

i dA 1 (r 1 _ g(*,r)M*,T) = 6(*,t)Q(*,t) 
c/r u; — w — 

CLT U — U0 2 US — UJ\ 

They take up stationary solutions 

A x {t) = A l (0) exp(-iA 1)2 r), A 2 (r) = A 2 (0) exj>(-i\ 1>2 r) (9) 
with characteristic values 

Ai = ^(p+VP 2 -4g), A 2 = ^(p- vV-V), (10) 

where 

p=l + + , q = • 11 

0J — LOi LO — UJ 2 UJ — L0\ 

These solutions should be normalized so that | (0) | 2 H- |A 2 (0)| 2 = 1. Thus we obtained 
adiabatic solutions, mapping the superposition of three bare states into a dressed state. The 
solutions corresponding to the characteristic values Ai and A 2 have the form 

r , , . A 1 (0,\ l! 2)M^)eM-iEa 1 t/h) + A 2 (0,\ 1 ,2)M^)eM-iEa 2 t/h + iT)\ 

+6(0, A 1)2 )^(^) exp(-iE b t/h) J 

(12) 

where 

Ai(0, Ai 2) = ; = =, (13) 

(W - W2)(Ai, 2 -q)y/l + ^6l67(^ - ^ 2 ) 2 (A li2 - g) 2 ' 

A 2 (0,Ai, 2 )- ' 



\A + -^ 2 ) 2 (A ll2 -g) 2 
Both excited state amplitudes 6(0, Ai) and 6(0, A 2 ) are determined from Eq.( 



III. LIGHT PROPAGATION IN QUANTUM PREPARED MEDIUM. SIMPLE 
APPROXIMATE ANALYTIC SOLUTION 



Now we proceed to the discussion of the field propagation in the medium using wave 
equation 

where the quantity pd (&(t)\ a \^(t)) jdt stands for the density of the bias electric current 
in dielectric gas medium, p is the number of atoms per unit volume and the brackets mean 
quantum mechanical plus thermal state averaging. The form above in fact assumes that 
the gas density is low enough and the temperature is high enough for quantum mechanical 
collective effects to be irrelevant. We make the ansatz 

(*(t)| ~i |*(t)> = 1 (a*(t, Ai) + p*(t, A 2 )| ~i \cfi!(t, Ax) + /3¥(t, A 2 )) + 

l + exp(-^fe^) 

(15) 

(m(t, Ai) + p*(t, A 2 ) I ~i \cM(t, Ai) + 0* (t, A 2 )) , 



+ ex P (ir 

which assumes independent formation of superposition states from the thermally populated 
levels 1 and 2. The coefficients (a, (3) and (a, /3) are the probability amplitudes corresponding 
to each channel. 

Standard transformations in frame of slowly varying amplitudes approximation bring us 
to a reduced form of wave equation (114j) . which for a better insight of the field dynamics we 
write in the following symbolic form: 

Tz ~ \b) v {Zl t) = ^ |2) + Gc ^ |2) cos((A2 " Al)r) + Gs ^ |2) sin((A2 " Al)r)) 71 {Z} T) 

(16) 

The wave amplitude rj (z, r) is defined in a dimensionless, symmetric with respect to both 
1 — 3 and 2 — 3 transitions form 

2d l d 2 E (z,T) 

" (Z ' T)= (d 1 + d ^A a - (17) 
The right-hand side expression in the brackets conditions the response of medium to the 
coupling field and the nature of the field evolution. Note that there is nothing surprising 
in the fact that the response of medium in superposition state is time dependent, as is seen 



from Eq.(16), which leads to phase modulation and relevant phenomena. The "hidden" 
from the first glance and crucial point here is the existence of an imaginary part in the 
response function for this non-dissipative medium under consideration. Periodic in time 
imaginary part acknowledges the quantum superposition nature of atomic state prior to the 
atom-field interaction and results in sequential amplification and suppressing of the field 
intensity during the propagation. And this is the reason (or origin) of the birth of series 
of short or ultrashort optical pulses, which we term as quantum superposition pulse train 
(QS-PT) generator. 

We will not discuss in this paper the questions connected with phase modulation and 
proceed to a solution of the problem in terms of the field intensity I77 t) | 2 . The wave 
equation for it follows from Eq.(16): 



d_ 



2 2ixptjjd\d<i 
= 1- 



( 



UJ — 0U 2 00 — 0J\ 



X 



l+exp( 



+ 



l+exp( 



a/3* 



{A 2 (0, Ai) A 1 (0, A 2 ) - A 1 (0, Ai) A 2 (0, A 2 )) 



exp (-i (A 2 - Ai) (C + t)) 
^exp(i(A 2 -A 1 )(c + r)) 



x 



(18) 



1 (>T 



with ( = A z/c, t = A (t — z/c). Eq.(18) is the main result of this paper and the basis of 
later discussion. We should note that in general it is a nonlinear equation and nonlinearity 
enters through the intensity-dependence of characteristic values Ai and A 2 . As we can see 
they appear in the equation in two forms: A 2 (0, Ai) A± (0, A 2 ) — A\ (0, Ai) A 2 (0, A 2 ) and 
A 2 — Ai. Numerical calculations, see Fig.2, show that the first one of them, denoted [A 2 , A±\ 
with high accuracy may be treated as not depending on wave intensity. The difference 
A 2 — Ai may also be viewed as as intensity independent, but only when i] (z,t) « 1 and 
equals to unity (A 2 — Ai = 1). For this intensity range the problem becomes to a great 
extend exactly soluble and has the following general solution form : 

/ 2jrpujdid2 ( _\ 1_\ x \ 



\v{z,t)\ 2 = |?7(0,t)| 2 exp 



(A 2 (0, Ai) A, (0, A 2 ) - A, (0, Ai) A 2 (0, A 2 )) x 

a*f3 1 a* /3 



1+exp 



l-cxp(-i(A 2 -Ai)(Ao2/c)) 



+ 



!+ ex p(^) 



X 



(19) 



v A2 _ Al ■ exp (-i (A 2 - Ai) Ao0) + c.c. J 

where the characteristic values Ai and A 2 are determined for the given initial wave-front. 
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FIG. 2: Functional form [^2,^1] (see the text) as function of scaled field intensity rj . Atomic 
transition frequencies lu% = 287351Ao and U2 = 287350Ao imply sodium atom's spectrum (see 
paragraph III of the main text), and input laser frequency oj = 287360Ao is detuned from these an 
order of magnitude larger than the doublet splitting Ao = oj\ — 1^2 ■ 

This formula shows all the immense advantage of the scheme under consideration: the 
output intensity is modulated at the frequency A 2 — Ai, which is the distance between the 
ground doublet levels. This is a remarkable and a prominent performance index of the 
presented scheme, since a very wide range of high stability repetition rates can be attained 
for the pulse trains by simply choosing different atomic species with doublet (or more) 
ground level structures. 

The second prominent advantage of the presented scheme is the presence of a prefactor 
1/Aq instead of mf\\*\ i ~E? m ^ ne original version [34], where A is detuning of the intense 
pump field relative to the two-level optical transition. To neglect the spontaneous emission, 
one usually has to work in the range where detunings are A > lO^Hz, which is at least one 
order of magnitude larger than A for alkaline metal gases. At the same time in old version 
we are limited with Q/A = 0.1 upper value in the additional coefficient. So the present 
scheme possesses much more freedom for such parameter values as gas density, and perhaps 
can be implemented even on weak optical transitions. 

Another preferable side of A-type interaction scheme is the acquisition of superposition 
between the closely posed ground states, as it is not a difficult task and may be implemented 
by means of magnetic fields as well. In case of a preparative radiation field, one has a* = a, 
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(5* = —(3, a* = —a and (5* = —(5. The expression (19) becomes somewhat simplified and 
the relative intensity \rj (z, t)\ 2 / \r) (0,t)| 2 (RI) takes the ultimate form 

^ 4iripu)d 1 d 2 ( _V ]_\ f a/3 . off A x \ 

Relative intensity= exp SA ° V™ ™J V 1+6X p(-^) J 

^ A,(0M)Ai(0,y-ArtM)M°M) (sin((A2 _ Al) (r _ Aoz/c)) _ sin((A2 _ Al) r)) J 

(20) 

Expression in the exponent periodically oscillates as a function of time and coordinate 
around the zero value. Therefore, in principle two types of field envelope transformations 
are possible: propagation with a periodic amplification and propagation with a periodic sup- 
pression. The answer to the question which process will be realized for a given superposition 
state is determined through the values of uj — uji and u> — u 2 differences. 



IV. NUMERICAL ANALYSIS OF PULSE TRAIN FORMATION 

For quantitative analysis we discuss the case of 3<Si/2— transition of 23 Na. The 

hyperfine levels F — 1 and F = 2 of ground level 3Si/ 2 constitute the doublet structure 
discussed. As the hyperfine splitting in the excited state 3P\/2 is orders of magnitude 
smaller than the 1771.6 MHz splitting of 3Pi/ 2 , we discuss the former as a collection of 
degenerate sublevels of same energy. We are also allowed to neglect all thermal transitions 
between ground levels F — 1 and F = 2 during the interaction. 

The incident wave is taken circularly polarized and this polarization is preserved during 
propagation through the medium. Transition matrix elements now should be calculated 
between the magnetic sublevels, and this should include selection rules of the total moment 
F and its projection M on the quantization axis aligned with the propagation direction. 
This, as usually, is realized using formulas (see, for example [36]) 



(F', M'\ d + \F, M) = (('F, M\ cL |F', M))*, 
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FIG. 3: Splitting of monochromatic incident wave intensity into a series of ultrashort pulses in a 
quantum prepared three-level atomic medium. Here zAq/c = tt , p = 4 x 10 12 cm -3 , T = W~®K, 



separating a projection-dependent prefactor from the reduced matrix element (F'\ d\F). 
Here d± = d x ± id y are circular components of dipole moment operator. Finally it was 
convenient to use the formulae 



introducing the oscillator strength f(F — > F') instead of reduced matrix element (F'\ d \F). 

The expression for the light field amplitude is of Eq.(19) form. The depth of modulation 
increases with atom density and thereby decreases the duration of each individual pulse 
in the sequence. Pulse duration may be significantly controlled using resonance detuning 
too. Fig. 3 graphically illustrates the situation, plotting the envelope of the field intensity 
as a function of time for various frequency detunings. Though the parameters are chosen 
for illustrative purposes, they are realistic for ordinary laboratory conditions. As might 
have been expected, pulse repetition rate is irrelevant to resonance detuning and mimics the 
doublet splitting frequency. As to pulse duration, is around hundred picoseconds for the 



a = (3 = 0.5, j3 = a 



i\/0.75 . All the other parameters are the same as in Fig. 2. 
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chosen frequency region and concomitant with power amplification. One will have the same 
picture for spatial evolution too. 

V. CONCLUSIONS AND OUTLOOK 

We showed that the gas of quantum prepared A-type atoms is very efficient in splitting 
of monochromatic wave into a sequence of ultrashort pulses, enlarging the spectrum of 
pulse generation mechanisms suggested or implemented so far. We emphasize that the 
repetition rate of ultrashort pulses does not depend neither on the light intensity nor its 
detuning from the atomic resonance. The sole quantity determining the repetition rate is the 
frequency separation between the ground doublet levels between which the superposition has 
been preliminarily established. Duration of each pulse, vice versa, is rigorously dependent 
on coupling parameters and effectively controllable. These regularities accompanied by 
mild conditions for the experimental realization presents the QS-PT generator as a viable 
alternative to more conventional methods. 

Finally we would like to note, that since the QS-PT generator is inevitably grounded on 
the quantum superposition principle, the studied phenomenon of abrupt shape-reformation 
of the classical electromagnetic field envelope may be regarded as an example of quantum 
macroscopic effect. 
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